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Abstract: Atmospheric modelling of the components of the visually close binary systems Hip70973 
and Hip72479 was used to estimate the individual physical parameters of their components. The model 
atmospheres were constructed using a grid of Kurucz solar metalicity blanketed models, and used to 
compute a synthetic spectral energy distribution for each component separately, and hence for the 
combined system. The total observational spectral energy distributions of the systems were used as a 
reference for the comparison with the synthetic ones. We used the feedback modified parameters and 
iteration method to get the best fit between synthetic and observational spectral energy distributions. 
The physical parameters of the components of the system Hip70973 were derived as: — 5700±75 K, 
= 5400 ± 75 K, log g a = 4.50 ±0.05, log g b = 4.50 ±0.05, R a = 0.98 ± O.O7i? , R b = O.89±O.O7.R , 
and 7T = 26.25 ± 1.95 mas, with G4 & G9 spectral types. And those of the system Hip72479 as: 
T c a ff = 5400 ± 50 K, Tjfa = 5180 ± 50 K, log g a = 4.50 ± 0.05, log g b = 4.60 ± 0.05, R a = 0.89 ± O.O7i? , 
R b = 0.80 ± O.O77? , and tt = 23.59 ± 1.00 mas with G9 & Kl spectral types. 



Keywords: stars: physical parameters, binaries, visually close binary systems, atmospheres modelling, 
Hip70973, Hip72479 



1 Introduction 

The Hipparcos mission revealed that many previously 
know n single stars were actually b i nary or multiple sys - 
tems (|Shatskii fc Tokovininl ir998: Balega et alll2002ft . 
Most of these resolved systems are nearby stars that 
appear as a single star even with the largest ground- 
based telescopes except when we use hi gh resolution 
techniques like speckle inter ferometry fSH ljBalega et al.l 
l2002l:lTokovinin et al.ll201CH and ada ptive optics (AO) 
([Roberts! l201ll : iRoberts et all l2005h . These systems 
are known as visually close binary systems (VCBS). 

The study of binary systems plays an important 
role in determining several key stellar parameters, which 
is more complicated in the case of VCBS. Hundreds of 
binary systems with periods on the order of 10 years or 
less, are routinely observed with high resolution tech- 
niques. In spite of that, there is still a paucity of in- 
dividual physical parameters for the systems' compo- 
nents. So, spectrophotometry with atmospheric mod- 
elling is a complementary solution to this problem, by 
giving an accurate determination of the effective tem- 
perature, radius and luminosity for each component 
of a binary system. The method was successfully ap- 
plied to some binary system like A DS11061, Coul289 , 
Coul 291, Hipll352 and Hipll253 dAl-Wardatl l2002al . 
l2007l . 120091 : lAl-Wardat fc Wid van! 1 20091 1. 

The two binary systems Hip70973 and Hip72479 
are well known VCBS. So, they fulfil the requirements 
to be analyzed by the aforementioned method in or- 
der to get their complete physical parameters. Ta- 
ble[T]contains basic data of the systems from SIMBAD, 
NASA/IPAC and The G eneva-Copenhagen sur vey of 
the Solar neighborhood (|Nordstrom et all [2004'). Ta- 



ble [2] c ontains dat a from Hipparcos and Tycho Cata- 
logues (|ESAlll997l V 

The system Hip70973 was discovered by Rossiter 
(1938.51) with the 27 inch ( 0.69 m) telescope a t the 
Lamont-Hussey Observatory l|Docobo et al.l200"oT). Or- 
bits o f the sy s tem h a d been calcul ated by ICouteaul 
([19601 1 . iMorell (|l970l 1. iHeintzl l[l98ll 1 (two orbits; the 
first one with period 45.4 yr and dynamical parallax 

0."020, and the second o ne with period 22.4 yr a nd dy- 

namic al parallax 0."029),[S6dcrhiclm] ([19991 1 and lDocobo et all 
(|2000l ) (elements of this orbit are listed in Table [3]). 

The system Hip72479 (ADS9397) was discovered 
by Aitken in 1916.40 at the Lick Observatory. Its or- 
bits had been calcu l ated by Ivan den Bosl (|l954l . Il945l . 

ll964l 1. lEggenl l| 19651 . 1 19671 1 fdifferent orbits us ing pho- 

tomet rical parallax 0."026l. [Soderhiehnl (|l999h and lDocobo et al.1 
(2000) (elements of this orbit are listed in Table [3]). 

The estimated parameters will enhance our knowl- 
edge about stellar parameters in general, and conse- 
quently help in understanding the formation and evo- 
lution mechanisms of binary stellar systems. 



2 Atmospheric modelling 
2.1 Hip70973 

We adopted the magnitude difference between the two 
components Am = 0™56 as the average of all Am 
measurements under the speckle filters 550nm/40 & 
551nm/22 (see Table |4| as the closest filters to the vi- 
sual. This value was used as an input to the equation: 



^ = 2.512~ Am , 

.72 



(1) 
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Table 1: Basic data of the systems 





Hip70973 


Hip72479 


ref. 




RST4529 


A2983 




«2000 


14 /l 31 m 00?650 


14 h 49 m 13?621 


1 


<^2000 


-05°48'08."46 


+10°12'52."06 


1 


WDS 


14310-0548 


14492+1013 


1 


Tyc 


4996-131-1 


921-918-1 


1 


HD 


127352 


130669 


1 


Sp. Typ. 


G5 


K2V 


1 


E(B-V) 


0.0499 


0.0268 


2 


Ay 


m 156 


0^083 


2 


l0gT eff 


3.722 


3.705 


3 


[Fe/H] 


-0.14 


0.11 


3 


M v 


4 m 75 


5 m 32 


3 



^IMBAD 
2 NASA/IPAC|http://irsa~ 



i pac.ca ltech.edu 
iNordstrom et al.1 (|2004h . 



From these we calculated a preliminary individual 
m v for each component as: m va — 8!" 19 and m v b ~ 
8 m 75. 



Table 4: Magnitude difference between the com- 
ponents of the system Hip 70973, along with niters 
used to obtain the observations. 



Am 


filter (A/AA) 


ref. 


m 42±0.15 


V Hp ■ 550wn/40 


1 


0™52 


550nm/40 


2 


0™55 


698wn/39 


2 


0™60 


551nm/22 


3 


0™60 


657n?7j/5 


3 



^ESAI (I1997D. ^Horch et all (120081) 



^Tokovinin et al.1 (|2010h . 



Table 2: 
logues 



Data from Hipparcos and Tycho Cata- 





Hip70973 


Hip72479 




HD127352 


HD130669 


Vj(Hip) 


7!™68 


87*42 


Bf 


8™667± 0.014 


97*540 ± 0.020 


V T 


77*781 ±0.011 


87*534 ±0.014 


(B - V)j(Tyc) 


0™775 ± 0.003 


07*866 ± 0.007 


ir Hip (mas) 


26.04 ± 1.04 


24.21 ± 1.29 


iiTyc (mas) 


30.6 ±8.7 


20.8 ± 10.6 


v*m v (mas) 


24.31 ±0.89 


22.59 ± 1.23 



(2007) 



Table 3: Orbital elements of the systems 
(|Docobo et alJl2000h 



Hip 


70973 


72479 


WDS 


14310-0548 


14492+1013 


P (yr) 


22.98 ±0.30 


9.98 ±0.04 




1993.62 ±0,02 


1988.059 ±0.03 


e 


0.499 ±0.010 


0.491 ±0.001 


a (arcsec) 


0.243 ± 0.002 


0.127 ±0.001 


i (deg) 


49.1 ±2.0 


45.8 ±2.0 


O (deg) 


13.8 ±2.0 


142.3 ±2.0 


uj (deg) 


121.0 ±2.5 


156.8 ±3.0 


TTdyn (mas) 


23.2 


21.1 



Teriastron transit time (yr). 



along with the visual magnitude of the combined sys- 
tem m v = 77*68 from Table [2] as an input to the equa- 
tion: 



Table 5: Magnitude difference between the com- 
ponents of the system Hip72479, along with filter 
used to obtain the observations. 



A, 



filter (A/AA) ref. 



0I™09 ± 0.72 

17*00 

07*40 

0™70 

07*60 



V Hp : 550wn/40 1 

657nm/5 2 

551nm/22 2 

657nm/5 2 

551nm/22 2 



^ESAI ()1997D . ^Tokovinin etaD (|2010h . 



Using t he following main sequ ence relations and 
tables (e.g.. lLandll99l lGravll20"u5h : 

M v = m v + 5 - 5 log(d) - A, (3) 
log(R/R Q ) = O.51og(L/L ) - 21og(T/T Q ), (4) 
log 5 = log(M/M Q ) - 2Iog(ii/.Rg) + 4.43, (5) 

we calculated the preliminary input parameters (bolo- 
metric magnitudes, luminosities and effective temper- 
atures) of the indivi dual compo ne nts. We use d bolo- 
metric corrections of lLand (|l992f ) fc lGravl ||2005D . T Q = 
5777K and extinction (A v ) given in Table[T]by NASA/ 
IPAC. 

These calculated input parameters allow construc- 
tion of model atmospheres for each component using 
grids of Kurucz's 1994 blanketed models (ATLAS9), 
where we used solar abundance model atmospheres. 
Hence a spectral energy distribution for each compo- 
nent can be built. 

The total energy flux from a binary star is created 
from the net luminosity of the components a and 6 
located at a distance d from the Earth. So we can 
write: 



-2.51og(/i+/ 2 ). 



(2) 



Fx 

from which 



<r = m 



Rl + H b x 



(6) 
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F x = (R 2 a /d) 2 (H a x +H b x -(R b /R a ) 2 



(7) 



where HI and H b x are the fluxes from a unit surface of 
the corresponding component. Fa here represents the 
total SED of the system. 

Within the criteria of the best fit, which are the 
maximum values of the absolute flux, the shape of the 
continuum, and the profiles of the absorption lines, 
and starting with the preliminary calculated parame- 
ters, many attempts were made to achieve the best fit 
between the observed flux and the total computed one 
using the iteration method of different sets of param- 
eters. 

Using Hipparcos modified parallax (ji — 24.31 ± 
0.89mas. Table [2}, the best fit was achieved using the 
following set of parameters: 



T c Q ff = 5700 ± 75K, T 6 ff = 5400 ± 75K, 

log g a = 4.50 ± 0.05, log g b = 4.50 ± 0.05, 
R a = 1.05 ± O.O7i? and R b = 0.96 ± 0.07^©. 

But the va lues of the e st imated radii disagree with 
those given bv lGravl (|2005l ) . iLanel [|l992D and the R-L- 
T relation (equation [JJl for the main sequence stars. 
According to equation [7] this disagreement refers to 
a misestimation in the parallax of the system, which 
means that changing the parallax of the system affects 
the values of the components' radii. 

So, in order to reach reliable parameters for the 
system, we went the other way; i.e. we started with 
the ra dii which are compatible with the tables of lGravl 
(2005) and changed the parallax. Keeping in mind the 
values of the total observational Vj , Bt , Vt and Am as 
our goal in achieving the best fit between the synthetic 
and observational total absolute fluxes, we reached 
that using the following set of parameters (Fig. [1]): 

T c Q ft = 5700 ± 75K, T* B = 5400 ± 75K, 

log g a = 4.50 ± 0.05, log g b = 4.50 ± 0.05, 
R a = 0.98 ± 0.07.R© , R b = 0.89 ± 0.07^© , 



and 



d = 38.10 ± 3.04pc(?r = 26.25 ± 1.95mas.), 



Thus the luminosities of the components follow as: 
L a = 0.91 ± 0.08L©, and L b = 0.61 ± 0.05L©. These 
values represent adequately enough the parameters of 
the systems' components. 

Fig.[T]shows the best fit between the tot al synthetic 
SED an d the observational one taken from I Al-Wardatl 
(2002t]). Note that some of the strong lines and de- 
pressions, especially in the red part of the spectrum 
(around A6867A, A7200A, and A7605A), are H2O and 

O2 telluric lines and depressions. 

D epend ing on the tables of [Gray (2005) or using 

ILanel (|l992h Sp — T e g empirical relation, the spectral 
types of the system's components can be estimated as 
G4 & G9. 




5000 6000 7000 8000 

Wavelength (A) 



Figure 1: Dotted line: the total observational 
SED in the continuous spectrum of the system 
Hip70973. Solid lines: the total computed SED 
of the two components with d — 38.10 ± 3.04 pc 
(n = 26.25 ± 1.95 mas), the computed flux of the 
primary component with T c g = 5700 ± 75 K, log 
g = 4.50 ± 0.05, R = 0.98 ± 0.07i? Q , and the com- 
puted flux of the secondary component with T e g = 
5400±75K, log g = 4.50±0.05, R = 0.89±0.07i? Q . 



2.2 Hip72479 

For this system we adopted the magnitude difference 
between the two components Am = 0'. n 50 which is 
the average of all Am measurements under the speckle 
filters 551nm/22 (see Table [6]) as the closest filters 
to the visual. This value when used as an input to 
equation [T] along with the total visual magnitude of 
the system m v — 8™ 534 from Table [2] as an input to 
equation [21 results in the preliminary individual m„ 
for each of the component as: m va = 9!™07 and m vb — 
9 m 57. 

Following the same procedures explained in the 
previous section, and using Hipparcos modified paral- 
lax measurement as initial input (n = 22.59± 1.23, d — 
44.26pc) as a temporary input for the calculations, 
the best fit between the total synthetic SED an d the 
observational one taken from lAl-Warda t (2002b) was 
achieved using the following set of parameters (Fig. 



EJ: 



T c a ff = 5400 ± 50K, Teg = 5180 ± 50K, 

log g a = 4.50 ± 0.05, log g b = 4.60 ± 0.05, 

R a = 0.92 ± 0.07R e and R b = 0.84 ± 0.07J? Q . 

Again here, the valu es of the estimated radii d o not 
fit exactly those given bv lGravl (|2005l ). ILanel (|1992| ) and 
the R-L-T relation (equation [4]) for the main seq uence 
stars. So, we recalculated the absolute flux using iGravl 
(2005) radii as postulate values leaving the parallax 
subject to change. Following the guidance of the total 
observational Vj, Bt, Vt and Am, the best fit between 
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the synthetic and observational total absolute fluxes 
was achieved using the following set of parameters 

T c a ft = 5400 ± 50K, T c fl = 5180 ± 50K, 

log g a = 4.50 ± 0.05, log g b = 4.60 ± 0.05, 
R a = 0.89 ± 0.07J?© , R b = 0.80 ± O.O7,R , 

and 

d = 42.40 ± 1.85pc(7T = 23.59 ± l.OOmas.), 

Thus the luminosities follow as: L a — O.61±O.O5Z/0, 
and Lb = 0.41 ± 0.03Lq, with spectral types G9 & 
Kl for the primary and secondary components respec- 
tively. 



Table 6: Magnitude difference between the com- 
ponents of the system Hip72479, along with filter 
used to obtain the observations. 



Am 



filter (A/ A A) ref. 



T ." 09 ± 0.72 V Hp : 550nm/40 1 

l m 00 657nm/5 2 

0™40 551nm/22 2 

0? l 70 657nm/5 2 

m 60 551nm/22 2 

^ESAl (11997D . ^Tokovinin etaD (l2010l) . 



Entire flux (observational) 
- Entire flux (theoretical) 




3 Synthetic photometry 

In addition to the direct comparison, we can check 
the reliability of our method of estimating the physi- 
cal and geometrical parameters by comparing the ob- 
served magnitudes of the combined system from dif- 
ferent ground or space based telescopes with the syn- 
thetic ones. For that , we u sed the following relation 
(jMai'z Apellanizll2006l . 120071 ) : 



m p [F x>s {\)] = 



fP P (A)^(A)AdA 
J ^p(A)f\,r(A)AdA 



(8) 

to calculate the total and individual synthetic magni- 
tudes of the systems, where m p is the synthetic mag- 
nitude of the passband p, -P P (A) is the dimensionless 
sensitivity function of the passband p, F\ tS (\) is the 
synthetic SED of the object and F\, r (A) is the SED 
of the reference star (Vega). Zero points (ZP P ) from 
iMaiz Aocllanizl (|2007l ) (and references there in) were 
adopted. 

The results of the calculated magnitudes and colour 
indices of the combined system and individual compo- 
nents, in different photometrical systems, are shown in 
Tables [7] & [8] 



Table 7: Magnitudes and colour indices of the 
synthetic spectra of the system Hip70973. 



5000 6000 7000 

Wavelength (A) 



Sys. 


Fil. 


total 


comp. 


comp. 






cr = ±0.03 


a 


b 


Joh- 


U 


8.74 


9.17 


9.96 


Cou. 


B 


8.43 


8.92 


9.55 




V 


7.68 


8.20 


8.73 




R 


7.28 


7.82 


8.29 




U -B 


0.31 


0.25 


0.41 




B-V 


0.76 


0.72 


0.82 




V - R 


0.41 


0.38 


0.44 


Strom. 


u 


9.88 


10.31 


11.11 




V 


8.84 


9.30 


9.99 




b 


8.09 


8.59 


9.17 




y 


7.65 


8.17 


8.69 




U — V 


1.04 


1.01 


1.12 




v — b 


0.75 


0.71 


0.82 




b-y 


0.45 


0.43 


0.48 


Tycho 


Bt 


8.63 


9.10 


9.77 




Vt 


7.76 


8.28 


8.82 




Bt — Vt 


0.87 


0.83 


0.95 



Figure 2: Dotted line: the total observational 
SED in the continuous spectrum of the system 
Hip72479. Solid lines: the total computed SED 
of the two components with d = 42 .40 ± 1.85 pc. 
(it = 23.59 ± 1.00 mas), the computed flux of the 
primary component with T c ff = 5400 ± 50 K, log 
g = 4.50 ± 0.05, R = 0.89 ± 0.07^0, and the com- 
puted flux of the secondary component with T c g — 
5180±50K, log g = 4.60±0.05, R = O.8O±O.O7i? . 



4 Results and discussion 

Looking deeply at the achieved best fit between the 
total synthetic SED's and the observational ones (Figs. 
[T]&[2]), we see that there is a good overall coincidence 
in the maximum values of the absolute fluxes and the 
shape of the continuum except for the blue part of 
the spectrum around A4700A (Fig. [3} . A part of this 
disagreement is due to the lack of some opacities in the 
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Tabic 8: Magnitudes and colour indices of the 
synthetic spectra of the system Hip72479. 



Sys. 


Fil. 


total 


comp. 


comp. 






cr = ±0.03 


a 


b 


Joh- 


JJ 


9.74 


10.19 


10.91 


Cou. 


B 


9.28 


9.78 


10.36 




V 


8.43 


8.96 


9.47 




R 


7.97 


8.52 


8.98 




U-B 


0.46 


0.41 


0.56 




B-V 


0.85 


0.82 


0.89 




V-R 


0.46 


0.44 


0.49 


Strom. 


u 


10.90 


11.34 


12.08 




V 


9.74 


10.22 


10.86 




b 


8.88 


9.40 


9.93 




y 


8.39 


8.92 


9.37 




u — V 


1.15 


1.12 


1.23 




v — b 


0.86 


0.82 


0.92 




b — y 


0.49 


0.48 


0.51 


Tycho 


Bt 


9.51 


10.00 


10.60 




V t 


8.52 


9.05 


9.56 




Bt — Vt 


0.98 


0.95 


1.04 



synthetic SED's and to the difference in the resolution 
between the synthetic and observational spectra. 




Table 9: Comparison between the observational 
and synthetic magnitudes, colours and magnitude 
differences for both systems. 



Hip70973. 


Obs.* 


Synth, (this work) 


Vj 


7™ 68 


7™68±0.03 


Bt 


8™667 ± 0.014 


8 m 63 ± 0.03 


Vt 


7!™781 ± 0.011 


7? l 76 ± 0.03 


(B - V)j 


m 775 ± 0.003 


m 76±0.04 


Am 


m 56 t 


m 53±0.04 


Hip72479 


Vj 


8 m 42 


8™43 ± 0.03 


Bt 


9™ 540 ± 0.020 


9 m 51 ± 0.03 


Vt 


8™ 534 ± 0.014 


8 m 52±0.03 


(B-V) j 


0™866 ± 0.007 


m 85 ± 0.04 


Am 


OT^O 1 


m 51 ±0.04 



5000 5500 6000 

Wavelength (A) 



*See Tabled 

t Average value for the filters 550nm/40 & 551nm/22 
(Tables H& gj) . 



lax is not excluded at all, especially if we kno w that 
the va lues of the radii were highe r in the table s of lGravl 
(120051 ) comparing with tables of lLand (|l992ft . 

Fig. [3] shows the posi tions of the c o mpon ents on 
the evolutionary tracks of iGirardi et all (|2000T t. where 
the error bars in the figure include the effect of the 
parallax uncertainty. The ages of the systems can be 
established from the evolutionary tracks. 

It is clear from the parameters of the system's 
components and their positions on the evolutionary 
tracks that they are solar type main sequence stars, 
in the early stages of their life. Depending on the 
formation theories, fragmentation is a possible pro- 
cess for the f ormation of th e systems studied in this 
work. Where iBonn ell (1994) concludes that fragmen- 
tation of a rotating disk around an incipient central 
prot ostar is possible, as long as t here is continuing in- 
fall. IZinnecker fc Mathieul (|200lh pointed out that hi- 
erarchical fragmentation during rotational collapse has 
been invoked to produce binaries and multiple systems. 



Figure 3: An expanded look at the fit between 
the total synthetic SED's and the observational 
ones for the left part of the spectrum of the system 
Hip72479. 

A comparison between the synthetic magnitudes, 
colours and magnitude differences with the observa- 
tional ones (Table [5]) shows a very good consistency 
within the error values. This gives a good indication 
for the reliability of the estimated parameters of the 
individual components of the system, which are listed 
in Tables [10] &[TT1 

Concerning the accuracy of our estimated paral- 
laxes, it depends on the accuracy of the input radii as 
it is clear from equation [6]and [7] This means that the 
first solution which depends on Hipparcos new paral- 



Table 10: Parameters of the components of the 
system Hip70973. 



Component 


a. 


b 


Teff (K) 


5700 ± 75 


5400 ± 75 


Radius (R©) 


0.98 ±0.07 


0.89 ±0.07 


logs 


4.50 ±0.05 


4.50 ±0.05 


L(L ) 


0.91 ±0.08 


0.61 ± 0.05 


m: 


4.97 ±0.10 


5.77 ±0.10 


Mass, (M Q )* 


1.07 ±0.08 


0.94 ± 0.05 


Sp. Type* 


G4 


G9 


Parallax (mas) 


26.25 ± 1.95 


Age (Gy) 


2.7 ±0.3 



'depending on the tables of lGravl (200 
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5. Because of the high similarity of the two compo- 
nents within each system, fragmentation is pro- 
posed as the most likely process for the forma- 
tion and evolution of both systems. 



1.0 




















' Hip 70973 1M ©^ 

T Component a 
- V Component b 

Hip 72479 

► Component a 

> Component b 


Main Sequence Qg"-^ 





3.80 3.75 3.70 3.65 

L °9 r eff 

Figure 4: The systems' components on the evolu- 
tionar y tracks (solar abundance) of iGirardi et al.l 
(l200(]h . 

Table 11: Parameters of the components of the 
system Hip72479. 



Component 


a. 


b 


Teff (K) 


5400 ± 50 


5180 ± 50 


Radius (R©) 


0.89 ± 0.07 


0.80 ±0.07 




4.50 ±0.05 


4.60 ±0.05 


L(L ) 


0.61 ±0.05 


0.41 ± 0.03 


M* 


5.77 ±0.10 


6.25 ±0.10 


Mass, (M©)* 


0.94 ±0.05 


0.85 ±0.04 


Sp. Type* 


G9 


Kl 


Parallax (mas) 


23.59 ± 1.00 


Age (Gy) 


2.7 ± 


0.3 



* depending on the tables of iGravl 1,2005). 



5 Conclusions 

The analysis of the two VCBS, Hip 70973 and Hip72479, 
using atmospheric modelling results in the following 
main conclusions. 

1. The parameters of the systems' components were 
estimated depending on the best fit between the 
observational SED and synthetic ones built us- 
ing the atmospheric modelling of the individual 
components. 

2. New parallaxes of the systems were estimated 
from these stellar parameters. 

3. From the parameters of the systems' compo- 
nents and their positions on the evolutionary 
tracks, we showed that the components within 
each system are similar solar type main sequence 
stars (G4 & G9 for Hip70973 and G9 & Kl for 
Hip 72479). 

4. The total and individual U BVR Johnson-Cousins, 
uvby Stromgren and BV Tycho synthetic mag- 
nitudes and colours of the systems were calcu- 
lated. 
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